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The intrinsic spin Hall effect (SHE) originates from the topology of the Bloch bands in momentum 
space. The duality between real space and momentum space calls for a spin Hall effect induced from 
a real space topology in analogy to the topological Hall effect (THE) of skyrmions. We theoretically 
demonstrate the topological spin Hall effect (TSHE) in which a pure transverse spin current is 
generated from a skyrmion spin texture. 


Transverse spin accumulation in semiconductors due 
to extrinsic spin-orbit scattering was first predicted by 
Dyakonov and Perel mm- Strong spin-orbit coupling 
(SOC) of the disorder scatters different spins in oppo¬ 
site directions leading to a non-zero transverse spin cur¬ 
rent perpendicular to the charged current. Evidences of 
the predicted asymmetric scattering of different spins was 
later abserved in optical (3] and photovoltaic |J experi¬ 
ments. Hirsch named this phenomenon the ‘spin Hall ef¬ 
fect’ (SHE) and proposed that the chargeless transverse 
spin current can be transferred back to a Hall voltage 
using an inverse SHE measurement [5]. Later theoretical 
studies predicted an intrinsic contribution to the SHE in 
the presence of SOC due to the topological property of 
the Bloch states at the Fermi surface EHg. Direct obser¬ 
vations of the SHE have been experimentally achieved in 
semiconductors using Kerr rotation microscopy |101 ITT] . 

In magnetic materials due to SOC, extrinsic or intrin¬ 
sic mechanisms can lead to a non-linear contribution to 
classical Hall signal fT21I14| . The non-linearity which is 
proportional to the magnetization is a result of the trans¬ 
verse accumulation of itinerant majority spins resulting 
in the anomalous Hall effect (AHE) |T5]. Similar to the 
SHE, the AHE can result from an intrinsic or extrin¬ 
sic mechanism. The intrinsic contribution to the AHE 
is related to the Berry curvature within the Fermi sur¬ 
face, which is determined by the topological nature of the 
Bloch bands [HI HR] , 

The momentum-space topological origin of the intrin¬ 
sic AHE is the same to that of the intrinsic SHE. Sim¬ 
ilarly, the real-space topology of a magnetic system can 
also induce a Hall effect m- An electron hopping 
through magnetic sites with particular chiral textures ac¬ 
quires a Berry phase and thus experiences an emergent 
gauge field during transport ESI- The emergent gauge 
field generates a Hall voltage that does not originate 
from SOC, which is usually referred to as the ‘topological 
Hall effect’ (THE) [T5]. Recently, a skyrmion lattice, a 
topologically non-trivial chiral spin texture, has been ob¬ 
served in helical magnets with a Dzyaloshinskii-Moriya 
(DM) interaction |20H22| . These materials provide ro¬ 
bust samples where the THE has been detected, and the 
measured Hall signal is a signature of the skyrmion phase 


in many B20 magnetic compounds [231126) . 

In the adiabatic limit, each electron spin passing 
through a single skyrmion has its spin aligned with the 
direction of spatial magnetization of the skyrmion which 
generates an emergent gauge field of up to one flux quan¬ 
tum [251 . This flux quantum confined in the area of a sin¬ 
gle skyrmion gives a gigantic effective field, that makes 
the THE a possible detection method for skyrmions. 
Moreover, the direction of the local magnetic field gener¬ 
ated by this emergent gauge field is opposite for parallel 
and antiparallel spin, which deflects them in opposite di¬ 
rections. This might separate the spin current from the 
charge current, generating an unconventional topological 
spin Hall effect (TSHE) which does not originate from 
band topology. Motivated by these possibilities, in this 
letter we theoretically investigate the THE and the TSHE 
resulting from a single magnetic skyrmion. The TSHE 
phenomenon discovered here can be explained in terms of 
a general physical picture that would apply equally well 
to a skyrmion lattice. 

Due to the lack of periodicity, we apply the non¬ 
equilibrium Green’s function method (NEGF) to simu¬ 
late the coherent transport of itinerant spins traversing 
a single magnetic skyrmion m- The tight-binding elec¬ 
tron Hamiltonian we employ is, 

H e = - J H X C ^ Ci ■ S ^ t J2 + h ' C ') ’ W 

* (*>i) 

where oy is the spin of itinerant electrons, Jh is the 
Hunds’ rule coupling, t is the nearest neighbor hopping, 
and S,; is the local magnetization. It has been pre¬ 
viously discussed that the external magnetic field does 
not contribute much to the Hall effect, therefore we ne¬ 
glect its effect on the electron by taking the hopping 
parameter to be real m- Thus, the Hall signal ob¬ 
served in the following calculations is purely from the 
emergent gauge field of the skyrmion. The spin tex¬ 
ture {S.;} contains a single skyrmion located at the cen¬ 
ter of a 4-terminal cross bar (as shown in Fig. [lj. 
This texture is fully damped using the Landau-Lifshitz- 
Gilbert (LLG) equation with the magnetic Hamiltonian 
"s = - JE {i , j} s,• s j-DE {iJ) tij■ S, XS ; -V h.,S„ 

Here J is the nearest neighbor exchange coupling, D is 
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FIG. 1. (color online) The geometry of a 31 x 31 tight-binding 
cross bar. The arrows denote the in-plane component of the 
magnetization texture of a single skyrmion. The color plot 
demonstrates the S z component. The four terminals are num¬ 
bered clock-wise. 


the DM interaction and ho is the external magnetic field 
perpendicular to the cross-bar plane. For simplicity we 
choose D = J = h.Q. Periodic magnetic boundaries are 
applied at the terminals, while open magnetic bound¬ 
aries are used for other boundaries along the cross-bar, 
which gives large in-plane magnetization components at 
the edges. The skyrmion at the center is generated by 
manually creating a unity topological charge and then 
relaxing the spin texture until the magnetic energy is 
stable. Details of the magnetic dynamical simulations 
can be found in Ref. [ 251 . 


For the electron transport calculation, semi-infinite 
boundary conditions for electron states are applied to 
the four terminals of the cross bar. Each terminal is 
assumed to be a thermal bath of carriers with chem¬ 
ical potential /i rn . The semi-infinite electrodes are in¬ 
cluded by adding self-energy terms, E m = Gg^t, to the 
terminal blocks of H c , where g R is the surface Green’s 
function of terminal m. The retarded Green’s function 
of the device region bounded by the terminals is given 
by G r = [el — H e — S m ] _1 . In the linear response 

limit, the zero-temperature terminal currents, J m , are 
given by I m = ( e/h ) J2 n T m,n 5 y n - Sfx n denotes the chem¬ 
ical potential shift due to an applied bias in terminal n, 


(^Mn — Mn Cp). T m ,n — Tr 


(p r' 


(m n ) 


is the transmission coefficient between terminal m and 


where G(L = G£t, and r m = it S m - Sj, 


At 


steady state, the charge current is conserved such that 
Tm.m = - T ™n- A Symmetric bias is applied be¬ 

tween terminals 1 and 3 , Sy\ = —Sfj ,3 = 8/1 = O.IJh- 
Enforcing J 2 = I4 = 0 in the Hall effect measure¬ 
ment, the transverse Hall voltage can be solved as = 
5 fi (P — Q) / (P + Q) and < 5 /xi = Sy (R — S) / (R + S), 


where 


( P = T21T41 + X21T42 + T21X43 + X24T41 

I Q = T23X41 + T23T42 + T23T43 + T24T43 

| R = T42T21 + T21X41 + T23T41 + T24T41 

l S = T42T23 + T21T43 + T23T43 + T24T43 


Thus, the topological Hall angle can be evaluated as 
tan 0 T h = E h /E x = (/Lt 2 - A^VOi - M3)• 

Once Sfim and I m are obtained, then the total ter¬ 
minal spin current, = x,y,z ), is evaluated 

from I^ a = |Tr [er a I™ q ], where cr a = 1 (g) a a is 
the extended Pauli matrix and I" eq is the terminal 


current operator I“ q 


2nh 


SGZ si - 


G R 6^ - <5S“G£] , SG^ = E n G R , n T n , n G £ m <W, 
and 5 S™ = T m (e F )Sii m . The intensity of the TSHE 
is described by the spin Hall angle, a renormalized ratio 
between the transverse spin current and the longitudinal 
charged current 


$TSH 


2e\ &xy 
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where I13 = h — I3, and lf 2 a = lf z — I 2 Z - 

First, we study the THE and TSHE for the case of pure 
spin injection. By setting t = 0.2Jh, the tight-binding 
band-width is smaller than the spin splitting given by 
Jh- Therefore no matter where the the Fermi level lies, 
the electron injection does not mix different spins. The 
0th an d 0 tsh f° r different positions of e F are shown in 
Fig. E a) and (b), respectively. The corresponding sur¬ 
face density of states (DOS) that determines the type 
of current injection at terminal 1 is shown in Fig. [ 2 jc) . 
When the surface DOS is zero, the Fermi surface lies in 
the spin gap, and injection is absent, so that both 0 th 
and 0tsh are suppressed to zero. As e F passes through 
the bands, pure spin injection gives a Hall angle up to 
± 0.2 indicating the expected THE. The corresponding 
value of 0tsh is within ± 0 . 005 . At e F = ± Jh, both 0th 
and 0tsh change sign. 

The sign change of the Hall angles can be explained 
by the spin and carrier-type composition of the injec¬ 
tion from the ferromagnetic contacts. For each transport 
channel, a one-dimensional tight-binding chain gives a 
negative cosine electron band dispersion, which has a sign 
change of the effective mass at the band center. The ef¬ 
fective mass (m*) is positive at the bottom band-edge, 
and becomes negative at the top. When an up-spin elec¬ 
tron with positive m* is injected from terminal 1, it is 
scattered to the “right” due to the effect of the emergent 
gauge field generated by the skyrmion. This is denoted 
as scenario (I) in Fig. id). Alternately when m* < 0 , 
an up-spin electron injected from terminal 1 is equivalent 
to a down-spin hole injected from terminal 3 . Since the 
spin scattering due to the skyrmion is anti-symmetric, 
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FIG. 2. THE and TSHE for the case of pure spin injection 
(t = 0.2Jh). The (a) Hall angle 0th and the (b) spin Hall 
angle 0tsh are shown as a function of ef ■ The surface density 
of states at terminal 1 is shown in (c). The four scenarios of 
different carrier-type and spin compositions are illustrated in 
(d). 

the down-spin hole is scattered to its “left” as denoted by 
scenario (II). In a multi-channel scenario due to the trans¬ 
verse confinement, the tight-binding band splits into sev¬ 
eral sub-bands. Thus, the number of the electron bands 
and the hole bands crossing the Fermi level changes at 
different positions of ef- As ef moves from the bottom 
band-edge to the band-center, the number of electron 
bands crossing ef decreases, while the number of hole 
bands increases as depicted in Fig. [2jc) . Right at the 
band-center, the electrons and holes are equal, indicat¬ 
ing an equal contribution from both scenarios (I) and (II), 
which leads to a cancellation of both 0th and 0tsh ■ Fur¬ 
ther increasing ef, scenario (II) starts to dominate such 
that 0th and 0tsh change sign. Similar arguments can 
be applied to scenario (III) and (IV) for the down-spin 
case (see Fig. id)). 

Semiclassically, the relative strength of THE to the 
TSHE can understood as a cancellation of the transverse 
electric field due to charge accumulation at contacts (2) 
and (4) with the gauge field of the skyrmion. In all 
these pure-spin injection scenarios, the spin current is 
carried by charge which leads to a transverse accumu¬ 
lation of charge resulting in a Hall voltage and hence 
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FIG. 3. THE and TSHE in the case of mixed spin injection 
(t = 1.5Jh). (a) and (b) demonstrate the values of 0th and 
0tsh for different positions of tF- The red dashed lines corre¬ 
spond to the case where the central skyrmion is removed, (c) 
is a plot of the surface DOS at terminal 1. 

a THE. Since the transverse electric field cancels the 
Lorentz force given by the emergent gauge field of the 
skyrmion, a continuous spin current is suppressed at the 
steady state, making the TSHE insignificant. However, 
an order-of-magnitude increase in 0tsh can be achieved 
in the case of mixed spin injection which we discuss next. 

To simulate mixed spin injection, the hopping term 
is increased to t = 1.5Jh such that the injection band- 
widths of each spin are enlarged and overlap in some 
range of ef- The calculated values of the 0th and 0tsh 
are shown in Fig. [3] along with the corresponding results 
in the absence of a skyrmion for comparison. For energies 
in the range of —4 . 5 Jh < ep < —Jh and Jh < £f < 
4 . 5 Jh, 0th vanishes to ~ 0 , whereas 0tsh increases by 
approximately an order of magnitude compared to the 
case of pure-spin injection. Additionally, in the energy 
range — Jh < ef < Jh, the Hall angle corresponding to 
the THE 0p is finite and roughly same order as that in 
the case of pure-spin injection. 

To explain the presence of the TSHE, we again re¬ 
fer to the four scenarios shown in Fig. id). Within 
— 4 . 5 Jh < Ep < — Jh, the transport is dominated by sce¬ 
nario (I)+(111) as shown in Fig. ic). In this case, both 
the spin-up and spin-down electrons are injected from 
terminal 1. Due to the presence of a skyrmion there ex¬ 
ists a topological Hall signal which produces a transverse 
electrical field, Eth ■ At steady state, the zero-current 
condition at terminals 2 and 4 requires eU^H = — F t 
and e£TH = — F j. satisfied simultaneously. Due to the 
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chirality of the skyrmion, the emergent field experienced 
by the up spin is opposite to that experienced by the 
down spin, which generates opposite emergent Lorentz 
forces on the two types of spins (F j~= —F j,). Therefore, 
the zero-current condition in the transverse direction can¬ 
not be satisfied unless Eth = 0. Although imbalanced 
spin injection occurs due to the ferromagnetic electrodes, 
the THE must be suppressed in steady state as long as 
the transport is dominated by the same type of carrier. 
Since there is no electrostatic field to balance the emer¬ 
gent Lorentz force, a continuous chargeless spin current 
is established. Similar explanations [(II)+(IV)] can be 
applied for Jh < cf < 4 . 5 Jh- 

When the transport is dominated by two different 
types of carriers with the same spin, the TSHE is sup¬ 
pressed, and the THE voltage becomes finite. In our cal¬ 
culations, this occurs when ep is within [—Jh, Jh]; anc ^ 
the transport is dominated by the scenarios (II)+(III). In 
this case the down-spin electrons and holes are injected 
from terminals 1 and 3, respectively. The electrons and 
holes are scattered in opposite directions and then accu¬ 
mulate at terminals 2 and 4, respectively. Since the same 
spin is assigned to opposite charges, a non-zero Eth de¬ 
velops at terminals (2) and (4) resulting in a finite THE 
with a vanishing TSHE. 

To further demonstrate the differences between the 
THE and the TSHE, we show the vector map of the spin 
current density Jg_ (r) and the corresponding color map 
of the charge accumulation in Fig. [4] The spin texture 
and the terminal numbering are the same as in Fig. ^ 
For the THE case shown in Fig. 0 a), £f = —0.05Jh and 
the transport is dominated by scenario (II)+(III). There 
is a net drop in the transverse chemical potential between 
leads 2 and 4. The Jg, vectors circulate symmetrically on 
either side of the skyrmion, generating no significant to¬ 
tal transverse spin current. This corresponds to the case 
where 0 th ~ —0.2 and 0 tsh ~ 0. For the TSHE case 
shown in [Fig. 0b)], the transport is dominated by sce¬ 
nario (I)+(III). The equal-potential contour of Sfi (r) = 0 
cuts all the way across the vertical bar indicating little 
charge imbalance between leads 2 and 4. In transverse 
leads 2 and 4 there is a net spin current directed from 
lead 2 to lead 4 giving a negative 0tsh ~ —0.05. 

The TSHE discussed here is of similar magnitude as 
the SHE in broadly used Pt thin films [22]. However, the 
physical mechanism giving rise to the TSHE is funda¬ 
mentally different from the one leading to the spin Hall 
effect in strong spin orbit coupled systems. In such sys¬ 
tems, the spin Hall effect results from the topological 
property of the Bloch bands in momentum space. In 
contrast the TSHE results from the topological property 
of the skyrmion spin texture in real space. The real- 
space topology exerts opposite emergent Lorentz forces 
on different spins, which can induce the TSHE. 

We thank the helpful discussions with Prof. Jing Shi at 
UCR and Dr. K. M. Masum Habib at Univ. of Virginia. 



FIG. 4. (color online) Vector map of Jg, (arrow plot) and 
the effective chemical potential distribution (color map) for 
(a) the THE and (b) the TSHE. A longitudinal applied bias 
of 5[ii = —Sfi 3 = O.IJh is applied. For the THE (a), the 
spin current symmetrically circulates on either side of the the 
skyrmion resulting in no net transverse spin current. The 
electron and hole accumulation induces an imbalanced trans¬ 
verse potential distribution. For the TSHE (b), the transverse 
chemical potential distribution is symmetric, and a charge-less 
spin current is established in the transverse direction. 
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